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^ : Abstract 

^ \ Methods of analysis of nucleon-nucleus scattering data have progressed markedly over the past 

50 years, yet many analyses of scattering data still use prescriptions specified at various stages of 
that progress. But the assumptions made or inherent in those analyses are poorly justified and 
most are no longer necessary to facilitate evaluations. We investigate some of those assumptions 
and highlight deficiencies and the uncertainties they may make in data analyses. 



(N 






> 

p 

o 
o 



X 



PACS numbers: 



*Electronic address: amos@physics.uniinelb.edu.au 



'f Electronic address: |stevenka@ac.uj.za 



I. INTRODUCTION 

An ultimate aim of a nuclear reaction theory is to probe details of the structure of the 
nuclei involved. There are many and diverse reactions and observables one might use in that 
quest but herein we consider just cases of nucleon scattering (elastic and inelastic), focussing 
upon questions that need be raised regarding methods of analyses used. The questions raised 
apply to various extent to theories for the scattering and reactions involving other nuclear 
projectiles. 

First we deal with the problem of elastic scattering of nucleons. It is of great regret that 
many assume that elastic scattering data are not of interest experimentally or theoretically. 
But almost all theories of nucleon-nucleus reactions presume a knowledge of the elastic 
scattering data for the system. It is commonly thought that such can be approximated by 
using some global, phenomenological, local optical model potential. But in doing so, there 
are uncertainties that may vitiate using a reaction theory for data analysis to achieve the 
prime purpose of assessing the structure of the target. 

For many years analyses of experimental data from reactions initiated by nucleons, by 
protons mostly, have been of interest, in part as the incident proton most strongly interacts 
with neutrons in the target. Analyses of proton scattering data then complement those 
of electron scattering form factors which are more sensitive to the proton distributions in 
nuclei. If, as well, static moments and electromagnetic transition rates are known, the body 
of data could well assess putative details of low energy structures of nuclei. However, the 
methods of data analysis used need be of comparable surety for any conclusions drawn to 
have sufficient credibility. But one must bear in mind that electromagnetic transition rates, 
electron scattering form factors, and nucleon scattering data can be sensitive to disparate 
aspects of the nuclear wave functions so that a model of structure may be adequate for some 
data analyses and not others. Now while the theories used to analyze electromagnetic data 
are well established, with the size, if not specific influence, of corrections to results found 
known, that is not the case with most theories used to analyze nucleon-induced reactions. 

To date analyses of data from nucleon induced reactions can only be made using significant 
approximations to a basic theory of such scattering. In doing so, guidance from some physical 
properties and physics principles are needed to justify whatever model theory of the reaction 
is used. It is an evolutionary process and one can only assess approximations used in the past 
against results from another developed theory in which such approximations have not been 
made or are of lesser influence. To date, of course, all applicable theories of nuclear scattering 
involve approximations. Still, by treating the system more assiduously, it is possible to 
understand why conclusions drawn from past analyses need be treated cautiously and not 
taken as definitively as authors may have indicated. 

For analyses of nucleon-nucleus elastic scattering data, the first theoretical approach was 
phenomenological, dating back to Bethe [l|. The form chosen for the associated optical 
model potentials was taken to refiect the expected density distribution of nucleons in each 
target nucleus. Subsequently, optical potential was sought microscopically, folding a free 
nucleon-nucleon interaction with a density profile for the target. The most well known of 
those approaches was that of Kerman, McManus, and Thaler (KMT) J2i] . Later formulations 
took into account medium corrections to the interaction between the incident nucleon with 
those bound in the target as well as specifically including the nonlocality arising from ex- 
change amplitudes associated with treatment of the Pauli principle. Those (7-folding models 
are discussed in detail in a review |3|. 



All of the models identified in the above, and various sets with intermediate steps (such 
as the JLM approach [^) have been, and continue to be, used to analyse scattering data. 
What we present is not intended as a justification of the most recent applications. Rather we 
seek to point out the inadequacies of the models formed in distinct stages of development, 
why some should no longer be used, and why the conclusions drawn from the use of many 
be treated with more caution. 



II. THE OPTICAL MODEL POTENTIALS 

Most optical potentials used are of local type, whether they be defined phenomenologically 
(usually with Woods-Saxon form) or by folding schemes. Typical of the latter is the JLM 
method as recently used (5|, in which the direct (local) term is found by folding an NN 
interaction with the nucleon density profile of the target. The imaginary term is specified 
phenomenologically. Some also account approximately for the exchange process arising from 
the indistinguishability of the projectile with the nucleons bound in the target by using an 
additional local term. But there are problems with all of these approaches, and especially 
in using the relative motion wave functions generated from local optical potentials. 

A. On fitting elastic scattering data 

A primary problem, and a rather insidious one, is that justification of an optical potential 
is founded upon finding a fit to elastic scattering data: cross sections and spin observables. Of 
course, a viable optical potential in a single channel theory of elastic scattering should indeed 
provide such fits. But does a fit to such data imply a proper (the "correct" ) optical potential? 
Many assume that the answer to the question is yes. Few, it seems, fully appreciate the 
limited validity of such a justification. 

A fit to elastic scattering data only requires that the scattering theory produce a set of 
phase shifts that, when used in the standard sums of Legendre polynomials, give a good fit 
to data. In no way is that a proof of the uniqueness of the phase shifts, and even less of 
the interactions used in the Schrodinger equations producing them. Of course, a "proper" 
optical potential must lead to a good set of phase shifts, and as well a fit to measured elastic 
scattering data, but in specifying the interaction, one must use as many principles of physics 
as possible to have maximal credibility. Basically, all the usual phenomenological and part 
phenomenological models that end with purely local forms for a nucleon-nucleus optical po- 
tential fail in that respect to some extent. To get elastic scattering phase shifts only requires 
knowledge of the relative motion wave functions and their derivatives at extreme separation 
of the nucleon and nucleus. How those wave functions attained the required asymptotic 
properties from solutions of Schrodinger equations is not defined by those numbers alone. 

Given a set of phase shifts that typically are necessary to have an excellent fit to data, such 
as with the elastic scattering of 200 MeV protons from ^^C, phenomenological potentials, 
potentials from inverse scattering theories, and those microscopically formed (^f-folding) all 
exist from which quality fits to the same data have been found. Some of those fits are 
associated with a chi-squared per degree of freedom (x^/F) values ~ 1. 
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FIG. 1: (Color online) Comparisons of the differential cross sections (top) and analyzing powers 
(bottom) from defined sets of phase shifts with those deduced by using their inversion potentials. 



1. An example of equivalent optical potentials 

The conventional approach, the phenomenological optical model potential approach, to 
analyses of elastic scattering data is an example of inverse scattering. Specifically it is an 
example of numerical inverse scattering: the fit to data determines values of the parameters 
deemed best. More generally, the inverse scattering problem for fixed energy scattering 
resolves to the following: Given an S-function (phase shifts) at a particular energy, find local 
central and spin-orbit potentials that reproduce that S-function when used in Schrodinger 
equations. One fully quantal method to do that is the modified Newton-Sabatier scheme [gJ 
which finds potentials of the Bargmann class. That approach has been used with sets of phase 
shifts (the "data") obtained from a microscopic (f^-folding) model of the elastic scattering of 
protons at energies of 65, 100, 160, 200 and 250 MeV from ^^C The cross sections and spin 
measurables found therefrom |3| are in very good agreement with actual measured data. For 
comparison, those phase shifts have been used as input to the inversion scheme J6| to find 
a set of local (Schrodinger) potentials ^]. The (complex, central and spin-orbit) potentials 
were then used in Schrodinger equations to find if the same set of phase shifts would be 
regenerated. The cross sections and analyzing powers from the (^-folding model are shown 
by the symbols in Fig. [T] while those resulting from use of the determined (local) inversion 
potentials are shown by the various dashed lines. Clearly the reproductions are very good 
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FIG. 2: The potentials defined by using a modified Newton-Sabatier inversion scheme with the set 
of phase shifts for 65 to 250 MeV protons on ^^C. 



with but a minor mismatch for 250 MeV at large [Q > 70°) scattering angles. 

The potentials obtained from modified Newton-Sabatier inversion calculations are dis- 
played in Fig. |2j Therein the central real, central imaginary, spin-orbit real, and spin-orbit 
imaginary terms are displayed in panels labelled (a), (b), (c), and (d) respectively. The 
potentials for energies 65, 100, 160, 200, and 250 MeV are portrayed by the solid, dashed, 
long-dashed, dotted, and dot-dashed curves respectively. The energy dependence result- 
ing with these interactions clearly is severe. The central potentials become less refractive 
and more absorptive with increase in energy. The spin-orbit interactions vary even more 
markedly. These local potential forms determined by quantal inverse scattering theory are 
distinctly different from the Woods-Saxon interactions usually taken with the phenomeno- 
logical optical. But they fit the "data" (and from the origins, the actual measured data) very 
well. Indeed, a criterion was that the inversion process must give "data" fits with measure 
X^/degree of freedom of C(l). 

The actual cross section data [8|] for the elastic scattering of 200 MeV protons on ^^C 
are numerous and of high quality. The basic Newton-Sabatier method has been used |9|] 
to ascertain local complex potentials from them. Two distinct inversion potentials have 
been found. They are shown in the right panel of Fig. [Sj the real parts in the top and the 
imaginary ones in the bottom segments. The loop was completed by using both potentials 
in Schrodinger equations whose solutions gave phase shifts, and then the scattering cross 
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FIG. 3: (Color online) The data from the elastic scattering of 200 MeV protons on ^^C compared 
with cross sections determined from two inversion potentials (left panel). The two potentials are 
depicted, as functions of radius in the right panel. The real and imaginary parts are shown in the 
top and bottom segments respectively with their bands of confidence indicated. Details are given 
in the text. 



sections that are compared with the actual data in the left panel of Fig. [31 The evaluated 
cross sections are practically indistinguishable from the data reflecting that the two results 
fit the data with values of x^/degree of freedom 1.006 and 1.018. The two calculated cross 
sections differ only in the region of the largest scattering angle, and then but slightly. 

The two distinctive potentials were found from two sets of pole- zero pairs that represent 
the S'-function that one needs first for the chosen inverse scattering theory. Details are given 
in Ref. (9|. The approach allows specification of bands of confidence on those potentials. 
They are shown in Fig. [3] by the dashed curves surrounding the potentials. Those bands, 
determined from an error analysis |9|, are to be interpreted as follows: should another 
potential exist associated with a fit to the data with the same x^/F, then there is a 60% 
chance that it would lie within that band. Clearly then, as the two cases shown are so 
different, even such good data are insufficient to delineate a potential uniquely. To these 
one might add any finely tuned phenomenological model form. Thus a claim that a nucleon- 
nucleus interaction is the correct one, based solely upon its fit to elastic scattering data, is 
not valid. A "correct" potential must be based upon more physical information. Simply 
linking its form to expected matter and charge distributions of the nucleus is not sufficient 
in this regard, however. 



B. The effect of the Pauli principle 



A problem with most single channel (optical model) potentials is that they are local in 
form. Those defined fully, or in part, phenomenologically usually take parametric forms 
and values in accord with some notions regarding the matter and/or charge distributions of 



nucleons in the targets. The charge distribution of the targets may be well reflected by the 
form factor from the elastic scattering of electrons but one cannot assume that the matter 
distribution is then necessarily as well-described. This is most evident when describing 
scattering from or by neutron halo nuclei [10|. Certainly the elastic scattering of nucleons 
does not directly justify parameter forms for the reasons discussed above. 

One may hold faith in the neutron distributions given by large basis models of the ground 
state structure of any nucleus at least inside and near to the nuclear surface. Indeed, their use 
in folding A^A^-interactions to form an NA optical potential lead to results in good agreement 
with elastic scattering observables. But, in so doing, unlike the case of electron scattering, 
one must also take into account effects of the Pauli Principle. The prime effect of that 
principle in deflning NA optical potentials in the microscopic approach built essentially on 
KMT theory |2j], is to make the interaction highly nonlocal. That nonlocality is associated 
with the exchange amplitudes arising from the indistinguishability of the projectile with 
nucleons in the target. The contributions from scattering in which the detected (emergent) 
nucleon originally was bound within the target are extremely important. They must be 
treated carefully, since the associated (exchange) amplitudes have a different momentum 
transfer proflle to the direct scattering ones. Even more critical is that, often, these exchange 
amplitudes destructively interfere with the direct ones. This has been demonstrated for a 
range of energies [ll| so only an illustrative example is included herein. Cross sections from 
the elastic scattering of protons from ^^C at energies of 185 and 200 MeV are shown in 
Fig. m Masking of such effects in the fltting procedures of local phenomenological forms is 
insufficient. As discussed previously, such potentials are necessarily not unique. 

The considerable effects of the exchange amplitudes indicate also that theories predicated 
upon just the charge and/or matter densities have a problem. The exchange amplitudes 
require input of the one-body density matrices which include, but contain more information 
than, the simple density. One needs explicit single nucleon wave functions in fact. The 
severe destructive interference noted also means that seeking equivalent local potentials to 
approximate such may be quite problematic. 



C. Medium effects with the basic NN interaction 

In the 1980's there was a watershed in the studies of NN and NA scattering. First, 
reliable A'^A^ scattering amplitudes were defined and detailed phase shift analyses |l2| made. 
at least for energies up to the pion threshold. Second the Nijmegen [l3|, Paris [ij], and 



Bonn [15| A^A^ potentials were developed to fit those A^A^ amplitudes and phase shifts. 
Also the status of microscopic NA optical potential theories was reviewed at a seminal 
topical workshop in Hamburg p^j. Finally, experimental studies in those years produced 



many and varied high quality A^A^ data sets for energies to 1 GeV [17|, adding impetus 
to implementations of the theories. Since that time reviews |3|, ll8| of the optical model 
for scattering have been made given the central role it plays in studies of nucleon-induced 



reactions 1191-21 



Standard now is the optical potential built upon non-relativistic multiple scattering the- 
ory with the A^A^ scattering amplitudes modified from the free A^A^ values. The result is an 
effective NN interaction. Those modifications are caused by the two nucleons interacting 
within the nuclear medium and are due to Pauli blocking and mean field effects for both 
projectile and bound state nucleons. In addition, there are other effects due to the convolu- 
tion of the A^A^ scattering amplitudes with target structure that require off-the-energy-shell 
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FIG. 4: (Color online) Differential cross sections (top) and analyzing powers (bottom) from the 
elastic scattering of 185 (left) and 200 MeV (right) protons with ^^C. The solid curves are the 
results found using a complete (^-folding model optical potential while the dashed curves are the 
results obtained by omitting the nonlocal exchange parts of those interactions. 



scattering amplitudes. Also of importance is the complete antisymmetrization (at least at 
the two nucleon level) of the A + 1 nucleon scattering system which leads to direct and 
knock-out (exchange) amplitudes for NA scattering. As illustrated before, the effect of 
such NA exchange amplitudes is not small at any energy and those amplitudes are a major 



source of nonlocality. Other medium effects [22| also can have some significance. Details 



of the forming the effective interactions from a chosen free A^A^ interaction are given in the 
review |3|. We have used both the Paris and the BonnB interactions in our evaluations of 
A^A^ t-matrices and A^A^ ^f-matrices; the choice gives only small differences in results. 

In Fig. O the effects that medium dependence of the A^A^ effective interaction have in 
the cross sections and analyzing powers are shown. Results for the elastic scattering of 65 
MeV protons from ^°C are shown in the left panel while those from ^osp^^ ^^^ given in the 
right panel. The data J23| are shown by the filled circles while the calculated results found 
using gf-folding potentials are depicted by the solid curves. The results shown by the dashed 
curves are from optical potentials having no density dependence in the effective interaction, 
i.e. formed by folding the free A^A^ t-matrices. The density matrices for ^"Ca used in 



forming the folding potentials were taken from SHF (Skyrme-Hartree-Fock) calculations [24 
while those for 208p]^ were from a simple packed shell model with oscillator wave functions 
{hw = 6.7 MeV). 

The role of density dependence in the A^A^ effective interactions is very evident, most 
notably in the analyzing power results. Density dependence clearly improves agreement 
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FIG. 5: (Color online) Elastic scattering cross sections and analyzing powers for 65 MeV protons 
scattering from '^'^Ca (left) and from 208p|-, (^i-jgi^t). The solid curves are the results of using the 
density dependent, NN (^-matrices in forming the optical potentials, while the dashed curves are 
the results found using the free NN t- matrices in the folding process. 



with the data in both cases. However, it is apparent that the ^osp^^ results are not in as 
good agreement with the data as are the "^^Ca ones. The issue is one of structure as is 
shown next in Fig. O In this figure, cross sections from the elastic scattering of 65 MeV 
nucleons from ^"^Pb are shown. The data from proton scattering [23] and from neutron 
scattering 25| are compared with results from gf-folding potentials found using two different 
models of the nuclear wave functions. The dashed curves are results given by a packed shell 
model with proton oscillators oi Jtjuj = 6.7 MeV as before, but hw = 7.25 MeV was used 
for the neutrons. This gave a neutron skin thickness (5* = 0.16 fm), comparable to that 
of an SHF calculation [24]. The solid curves in the figure result when SHF model wave 
functions are used to specify the (/-folding potentials. The neutron scattering cross sections 
are quite well reproduced by both calculated results but the proton scattering results indicate 
a preference for the SHF model input. This reflects primarily the difference in the neutron 
distributions given by the two models of structure. 

Of course, by judicious adjustment of parameter values in a purely phenomenological 
approach, even better fits to this data can be achieved. Such is shown in Fig. [71 The 
parameter values used were those defined by van Oers et al. 26| for 40 MeV protons incident 



on ^°^Pb, and by Finlay et al. [27] for 40 MeV neutron scattering. Results found using the 
van Oers et al. potential for both scatterings are depicted by the solid curves while those 
found using the Finlay et al. potential are displayed by the dashed curves. The parameter 
values differ significantly as one may expect for nucleon scattering from N ^ Z systems. 
For completeness those parameter values are given in Table [B 
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FIG. 6: (Color online) Cross sections for the elastic scattering of protons (left) and of neutrons 
(right) from ^^^^Pb. 

TABLE I: The parameter values for the optical model potentials that lead to the fits to cross-section 
data shown in Fig. [71 Other specifics are given in Refs. [26] and [27] 



term 



40 MeV protons 



40 MeV neutrons 



Vq (MeV) 
ro (fm) 
oo (fm) 



51.8 
1.159 

0.784 



36.73 

1.205 
0.685 



Wq (MeV) 
Wd (MeV) 
ro (fm) 



3.98 
5.75 
1.321 



5.2 

2.55 

1.283 



aj) (fm) 


0.727 


0.569 


Vs.o. (MeV) 


6.66 


5.75 


Ts.o. (fm) 


1.044 


1.105 


as.o. (fm) 


0.905 


0.499 


Tc (fm) 


1.05 





D. Energy dependence problems 

The projectile energy dictates what approximations may be made to a 'full' scattering 
theory. At low values of incident energies, discrete state effects are known to influence 
scattering J28|. In those cases, a coupled-channel model of scattering is essential. An appro- 
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FIG. 7: (Color online) Cross sections for the elastic scattering of 40 MeV protons (left) and of 40 
MeV neutrons (right) from ^''^Pb, compared to results from using phenomenological optical model 
potentials. 



priate one, which ensures that the Pauli principle is not violated even with a collective model 
prescription for the coupling, now exists and has been used to explain compound nucleus 
structure even in exotic, radioactive light mass systems J29|. Treating the Pauli principle 
adequately in coupled-channel calculations was crucial as shown in a comparative study of 
two coupled-channel evaluations 



30 



Consider now, incident energies above any excitation energy value of giant resonances 
(if such exist in the target). Some time ago, it was shown that specific properties (the 
giant resonances) could dynamically influence proton scattering for energies of protons that 
coincide with their virtual excitation in the target [31.]. For such incident energies there 



has been another coupled-channels method proposed 32] • That method, the continuum 
discretised coupled channel (CDCC) method, has been used with some success in analyses 
of break-up reactions for projectile energies similar to those we have used. Those successes 
might be considered indicative of a problem with the ^f-folding approach for nucleon elastic 
scattering. However, the discretization of the continuum is only a statistical one and is not 
based upon any attributes of the target continuum. Also the effects of the Pauli principle 
(projectile-bound nucleons) are not taken into account with the phenomenological optical 
potentials used. But there are other concerns with the CDCC method. As Deltuva et al. |33| 
noted: "CDCC results obtained with two different basis sets, namely the basis set using the 
continuum of the projectile in the entrance channel, and the one using the continuum of the 
composite system in the final transfer channel, led to substantially different breakup cross 
sections for p(^^Be,^^ Be)pn. These findings raise concern about the accuracy of the CDCC 



11 



method, not only as a means to describe reaction dynamics, but also as an accurate tool to 
extract structure information on halo nuclei. " 

For medium energies, such as we consider herein, then such continuum couphng as well as 
coupling between the low excitation energies in the target are not expected to be important. 
Certainly there has been no need for such when a good model of structure and a reasonable 
A^A^ force were used in evaluations. That is so at least for elastic scattering cross sections 
usually greater than about 0.1 mb/sr. We do not dispute a role of coupled channels in a 
scattering process, but we are convinced that such are a requirement in analyses of nucleon 
elastic scattering data only when there are specific (collective and not too spread) states in 
the target nucleus at excitation in the vicinity of the incident energy value. Likewise there 
can be other reaction processes contributing to a given transition, such as the formation 



of a virtual deuteron via a p-d-p chain as considered by some [3J, |35j. But such second 
order terms, even if it is possible to effect a virtual deuteron and given the multiple final, 
possibly unbound, states to which the break-up could lead, should give small contribution 
in comparison to the direct ones in forming elastic scattering cross sections. Further the 
other problems as detailed before with either the modified JLM model or the CDCC coupling 
scheme that has been used 3J,|36| make the estimates given for the virtual formation process 
unreliable. These problems detract from conclusions made in those papers on the structures 
of the nuclei involved. 

Finally, as the incident nucleon energy increases to near and above the Delta resonance 
threshold in the NN frame, both that and other resonances need be taken into account in 
defining an effective A^A^ interaction. Relativity also becomes an increasing issue, and one of 
more than a kinematic factor. However, below about 300 MeV, the non-relativistic ^f-folding 
results compare well with data and more than favorably with global Dirac potential model 



results 37 



III. INELASTIC SCATTERING OF NUCLEONS 

Inelastic nucleon scattering data from excitation of discrete states in nuclei usually has 
been analysed using a DWA (distorted wave approximation). Coupled channels analyses also 
have been made, but they still are made using collective models of structure, usually with 
local coupled channel potentials. That localisation suffers from the same problems as does 
the local phenomenological or semi-phenomenological optical model potentials described 
previously. 

For elastic scattering one only requires appropriate relative motion wave functions at large 
distances; appropriate meaning that they yield a set of suitable phase shifts. For non-elastic 
reaction evaluations, with the DWA, the relative motion wave functions need be specified 
through the nuclear volume. The choice of optical model form then is of serious concern. 

A. Relative motion w^ave functions 

The problem of not considering the effects outlined in the previous sections (Pauli block- 
ing, mean-field effects, and nonlocality) extends to descriptions of nucleon-induced reactions 
for which the optical potential serves to define the effective interaction mediating the re- 
action. In that respect, an appropriate (physically justified) determination of the relative 
motion wave function (nucleon incident on the nucleus) is required in order to have confi- 
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dence in the description of the related reaction. In this subsection, we consider the relative 
motion wave functions associated with two optical model potentials. Specifically we com- 
pare the partial wave functions generated from a ^f-folding model [28|] and those from a 
phenomenological potential whose parameters have been set to 'best fit' the same cross sec- 
tion data. As shown above in Fig. El (^-folding model evaluations gave quite good fits to 
both the proton and neutron scattering cross sections. Even better fits to the data were 
obtained using phenomenological optical model potentials (see Fig [7]). Essentially then both 
the microscopic and phenomenological potentials yield relative motion wave functions that 
are the same asymptotically. However, their relative motion wave functions differ through 
the volume of the nucleus. 

Whatever be the (non-relativistic, spherical) nucleon-nucleus optical model potential, 
with the beam direction as the z axis, relative motion wave functions for neutrons are 
written in partial wave form 

vl/,(k, ^) = \Y. ^4^(2^ + 1) ^'A,(fcr) A 1 z. I J z/\ 3^r./a,), (1) 

where y is the standard tensor spherical harmonic for a spin-| function. 

Consider the case of 40 MeV neutrons elastically scattered from ^"^Pb. The relative 
motion wave functions for the s-, p-, and d- partial waves are displayed in Figs. El HI 
and [To] respectively. The real and imaginary parts of those wave functions are shown on 
the left and right hand sides of the figures. The solid curves portray the wave functions 
generated by the nonlocal, microscopic (yf-folding optical potential while those found using 
the phenomenological potential [27| are shown by the dashed curves. All wave functions were 



obtained from specific runs of Raynal's updated version, DWBA98, of his code DWBA91 |38 

The effect of nonlocality upon s-wave relative motion wave functions through the nuclear 
volume, shown in Fig. |H1 is to markedly decrease the real part while having a smaller decrease 
on the imaginary part of the wave function compared to that found using the (essentially 
phase equivalent) local potential model. The scaling is not the overall 10-15% that has often 



been assumed (the Percy effect [39|, |40|). 

The p- and d- wave functions are similarly effected, though in the case of the p-waves, 
it is with the imaginary part as given by the code, that the large size reduction due to 
nonlocality occurs. 

To emphasise the difference between the wave functions associated with the two 'phase 
equivalent' interactions for 40 MeV neutrons on ^'^^Pb, the partial wave sum completed along 
the beam line (for which fi^ = vr for z < and thereafter, gave the magnitude of the spin 
up {v = |) scattering waves as displayed in Fig. [TTl 

As nuclear structure models, by virtue of the fact that they are models, do not precisely 
describe nuclei, there are mismatches between observed quantities and the values predicted. 
Whether that be 7-decay rates or inelastic scattering cross sections, a match is sought by 
using effective transition operators with the ensuing scaling of calculated results defined as 
'core polarisation', which are considered in the next section. Nevertheless, the differences 
obtained in the relative motion wave functions when one considers the effects in scatter- 
ing (Pauli blocking, nonlocality, etc.) illustrate the problem in describing nucleon-induced 
reactions and the correction factors that may be necessary a posteriori. Such corrections 
were postulated as necessary in a recent article on coupled-channel descriptions of reaction 
cross sections of excited states which ignores the effects of Pauli and nonlocality in the 



specification of the effective interaction [41 . 
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FIG. 8: (Color online) s-wave relative motion wave functions for 40 MeV neutrons on ^^'^Pb. The 
real and imaginary parts are as indicated with the solid and dashed lines giving the results found 
using the ^-folding and local phenomenological optical model potentials respectively. 
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FIG. 9: (Color online) p-wave relative motion wave functions for 40 MeV neutrons on ^'^^Pb. The 
results for spins ^ ^-iid | are in the top and bottom segments respectively. Other specifics are as 
for Fig. El 

B. Core polarisation 

Core polarisation enhancement of calculated cross sections, or form factors, to replicate 
data is a measure of inadequacy of the structure model used. Such enhancement should 
be momentum-transfer dependent. That makes using any specified core polarisation from 
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FIG. 10: (Color online) d-wave relative motion wave functions for 40 MeV neutrons on ^'^^Pb. The 
results for 
for Fig. El 



results for spins | and | are in the top and bottom segments respectively. Other specifics are as 



one data set not necessarily the correct one for use with another. For example, to use 
electromagnetic transition rates (essentially zero momentum quantities) to specify a core 
polarisation for scattering cross sections typically specified for a range of momentum transfer 
(~ 0.2 to ~ 2.5 fm~^) is erroneous, as is the reverse. 



1. Core polarization from electromagnetic properties 

Consider the transition between the ground and first excited states in ^Li for which the 
B{E2) and longitudinal electron scattering form factor are both well known. However, first 
we need specify in brief the shell models to be used in evaluations of both quantities. More 



details are given in Ref. [42 



Consider shell model wave functions within the 0^, (0 + 2)hu, and (0-1-2-1- A)hu model 
spaces. The choice of model space dictates the choice of interaction. The ones used were 
the Cohen and Kurath (6-16)2BME interaction (CK) for the complete Ohw model space, the 
MK3W interaction for the complete (0-1-2) ^ shell model space, and the G- matrix interaction 
of Zheng et al. 43|] for the complete (0+2-1-4)^ model space. The shell model wave functions 
found with each set were used to find the longitudinal inelastic electron scattering form factor 
to the 3"*"; (2.185 MeV) state, amongst others |42|. In Fig. [T21 the results found using the 
(0 + 2 + 4)/ia; wave functions are depicted by the sohd curves, those with the (0 + 2) /to; 
ones by the dashed curves, and those with the Qhw model by the dot-dashed curves. In the 
panel (a), the data are those of Bergstrom et al. [46!] (circles), of Yen et al. [45|] (squares), of 
Bergstrom and Tomusiak 4^ (crosses), and Hutcheon and of Caplan 47| (triangles). The 
form factor, calculated with all shell models, is dominated by the C2 component, while the 
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FIG. 11: (Color online) The magnitudes of the complete scattering waves for 40 MeV neutrons on 
^'^^Pb along the z axis. The (^-folding interaction gave the result shown by the solid curve while 
that found from the phenomenological optical model is depicted by the dashed curve. 

C4 component is found to be negligible. With the MK3W and Zheng models, the calculated 
results reproduce the magnitude of the measured form factor above 1 fm~^. Both have 
strength from transitions outside of the Op shell which enhance the C2 strength. Such are 
missing in the Ofio; model. The (0 + 2 + 4)^ model structure is most favoured as there is 
almost exact agreement with the data in that region of momentum transfer. However, the 
B{E2) value for the 7-decay of this 3^; state is 9.3 ±2.1 e^fm^, and the values obtained by 



calculation using the three models of structure are significantly smaller [42j. So far as the 
7-decay is concerned, all calculations require a substantial renormalization to reproduce the 
measured value. That is confirmed by predictions of the electron scattering form factor at 
low momentum transfer. Below 1 fm~^ all of the calculated results are less than observation. 
Yet that degree of renormalization is not suggested by the results of the calculations of the 
form factor at higher momentum transfer. While this suggests that the internal (nucleon) 
dynamics of the nucleus are well described by the inclusion of higher hw excitations in the 
model space, such cannot account for the asymptotics of the structure. At large radii, which 
most influence scattering at low momentum transfer, the clustering of the wave functions 
are not reproduced by the shell model in which up to Ahw excitations are included. That is 
also why the shell model wave functions also do not give the correct quadrupole moment for 
^Li. This deviation of all the calculated results away from the data is illustrated further in 
the (b) panel of Fig. [H in which the B{E2 I, q) value for the 3+; (2.186 MeV) state in ^Li 
is shown as a function of momentum transfer. These values have been determined from the 
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FIG. 12: (a) Longitudinal inelastic electron scattering form factor to the 3~'';0 (2.186 MeV) state 



in ^Li, and (b) the B{E2 \, q) values, in units of e^fm^, as obtained from the form factor 44]. The 
data of Bergstrom et al. |44] (circles), of Yen et al. [45l] (squares), of Bergstrom and Tomusiak 
4g] (crosses), and of Hutcheon and Caplan [47] (triangles) are compared with the results of the 
calculations. The B{E2 X) value from the associated 7-decay rate is displayed by the diamond 
data point in (b). 



measured and predicted longitudinal inelastic scattering form factors, achieved by removing 
from the form factor most of the dependence on the momentum transfer according to the 
transformation of Brown, Radhi, and Wildenthal 48|. The B[E2 l) value is related to the 
associated 7-decay which is given by the g ~ intercept. 



2. Core polarization from radial transition form factors 

Core polarization has been required in DWIA (distorted wave impulse approximation) 
analyses of pion inelastic scattering to match calculated cross sections to observed data J49|. 
Therein, large basis microscopic models of nuclear structure were used to specify the 'col- 
lective' form factors for inelastic scattering to the 2^ states in ^^C and ^*Si. The 7r-nucleon 
t-matrices input were fixed by fits to vr — A^ elastic scattering phase shifts. All details are 



given in Ref. [49 



The large basis models of nuclear structure used in Ref. [49[, though superseded now, 
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FIG. 13: (Color online) Transition radial densities for the excitation of the 2^ {AAA MeV) state in 
12c (left) and of the 2+ (MeV) state in ^SSi. 

suffice for the point to be made. For the transition in ^^C, a large basis projected Hartree- 
Fock calculation was made. For the transition in ^^Si, large basis Nilsson model wave 
functions were used. All single nucleon orbits to Ahu excitation were considered and, from 
those model structures, transition one-body density matrices (OBDME) were defined. Of 
those the dominant set of values were in the p-shell for ^^C and in the sd-shell for ^^Si. 
Those values are very similar to what was found using the conventional Ohw shell model for 
each nucleus. However, the other entries were not trivial and can be considered as the 'core 
polarization' corrections to the basic Ohu model values. 

With those OBDME, and with harmonic oscillator single nucleon wave functions, radial 
transition densities for the excitation of the 2^ states from the ground states were deter- 
mined. They are the transition form factors required in the DWIA evaluations of the pion 
(inelastic) scattering cross sections. These radial transition densities for both nuclei are 
shown in Fig. [T31 

The solid curves are the transition densities found using the complete set of OBDME 
from the large basis wave model wave functions. The dashed curves are the densities formed 
by restricting the sums of OBDME to just the p-shell terms for ^^C and just the sd-shell 
entries for ^'^Si. Comparison of the two densities for ^^C indicates that the 'core polarization' 
enhances the density profile so that any cross section (essentially the Fourier transform of the 
transition density) would be simply increased in magnitude consistent with the conventional 
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usage of polarization charges. In fact the increase is a scale of ~ 1.4 so that a cross section 
would be increased by a factor ~ 2. However, that is not the case with ^^Si. Clearly the 
addition of the 'core polarization' contributions changes the profile of the transition density, 
reducing its strength in the region < r < 2 fni but increasing the strength at larger 
radii. The Fourier transforms of the two then will differ, and one can expect so also will the 
deduced cross sections. 



3. The Mn/Mp approach 



There are a number of studies, that of Chien and Khoa 50| being a most recent example, 
in which neutron scattering amplitude scalings were linked to proton ones through electro- 
magnetic properties. Those studies have used the Mn/Mp ratio defined by Bernstein, Brown, 
and Madsen 5l| who recognised that the combined measurements of inelastic hadron and 
of electromagnetic excitation of the same states in a nucleus could provide a way to probe, 
simultaneously, the proton and neutron excitations in the transition. They formulated a 
phenomenological link to specify the ratio of multipole transition moments for protons (Mp) 
and neutrons (M„), namely 



Mn 

M„ 






(2) 



where for the nucleus {Z, N), 6h is a deformation length to be determined from {p,p') scat- 
tering, 6eTn is the electromagnetically defined deformation length, and 6„, bp are interaction 
strengths of the incident proton with target neutrons and protons respectively. 
The multipole transition elements themselves are defined by 



M 



p(n) 



Jf 



E ^'^^(^^ 



iep{n) 



Ji 



pf''\r) r'^+Ur, 



(3) 



where Pxv') are transition densities as described in the preceding subsection and A is the 
multipolarity of the transition. 

Usually the B{E2) values are used 52|, |53|, though in Ref. |50|, it was the B{E2) values 
from the 2f state decays from proton-rich ^^Ne and ^°Mg to compare with their mirror 



systems in the neutron-rich Oxygen isotopes. Iwasa et al. [53| noted that there could be 
systematic uncertainties, possibly caused by insufficient knowledge of the mechanism of 
proton inelastic scattering. As the proton multipole transition moment is determined solely 
by electromagnetic data, the approach negates any contributions to Mp from the effective 
NN interaction. The process is not a reliable way to assess details of nuclear structure. 
Likewise it should not be used to specify B[E2) values from inelastic scattering cross sections 
from hydrogen targets, as has been considered for the spectra of radioactive ions. 

Though many of the problems with this approach have already discussed, it is worth- 
while listing them for emphasis. First the B{E2) are essentially a zero momentum transfer 
property which, if a perfect model of structure were available, would depend upon, only the 
proton transition density, with 



B{E2) 



i M^ 

{2J^ + l) P 



(4) 
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in units of e^fm^. It is questionable therefore to use this to compare data from a mirror 
system, for example the B{E2) from ^°Mg to that in ^°0. Of course, there are mirror 
symmetries in spectra of mirror nuclei, but save for the weak Coulomb displacement energies, 
the symmetry is due to charge independence of the nuclear force. However, at first order, 
the neutrons in a nucleus do not contribute to B{E2) values. They do contribute in so far 
as they influence proton distributions through the A^A^ interactions being stronger for a pn 
pair than with like nucleons, but not directly otherwise. 

Hadron inelastic scattering, of protons in particular, have measured cross sections usually 
in a finite range of momentum transfer values that are not very close to zero. That is 
particularly so with the inverse kinematic problem of radioactive ion beams (RIBs) scattering 
from hydrogen targets. As noted previously, there is a better match in momentum transfer 
values if one compares electron form factors with {p,p') cross sections, but experimental 
data from electrons scattering with RIBs, either a cross beam by or using a SCRIT (self 
contained radioactive ion target) experiment |5J] are yet to be realised. 

As stated earlier, perhaps the biggest problem of using the Mn/Mp ratio method is that 
it assumes that the proton scattering transition moment may be solely determined by fits 
to electromagnetic data. The allowance for specific nuclear interaction effects, the role of 
the Pauli principle in making the actual transition amplitudes creating exchange amplitudes 
which often have serious interference with the direct ones, and the effects of nuclear distortion 
in the relative motion wave functions, are all missing in the identifications given in Eqs. ([2]) 
and ([3]). These reasons bring into question the conclusions made in the recent article ^] 
dealing with neutron transition strengths of 2^ states in the neutron-rich Oxygen isotopes. 
That is especially so since the assessment of the structure of those isotopes specified by a 
large basis shell model and their use in analyses of the scattering cross sections of those 
RIBs from hydrogen targets existed in the literature ^] . 

IV. CONCLUSIONS 

Of course, if only a set of elastic scattering phase shifts that lead to a quality fit to elastic 
scattering and total reaction cross-section data are required, then any prescription that does 
so can be used. But usually the elastic scattering data and its fit are but a prelude to what 
is sought. By itself, one might seek to make conclusions about matter distributions in a 
nucleus from an assessment of elastic scattering data. To do so, however, one need use a 
more detailed theory of that scattering than usually has been the case. At least the Pauli 
principle and the nonlocality of the NA optical potential that it creates must be taken into 
account, and those nonlocalities should not be approximated with some local equivalent 
form. There is no need to do so now as the suite of programs DWBA98 and above, facilitate 
such analyses. The effect of medium corrections in the projectile nucleon - target nucleon 
effective interactions can be handled in those programs as well, and a utilitarian prescription 
for those effective interactions has been published [3|]. The computer programs that generate 
the effective interactions by those means are available upon request. 

Often the associated relative motion (distorted) wave functions are required, such as in 
DWA analyses of inelastic, charge and particle exchange reaction data. Those analyses seek 
conclusions regarding structures of excited states and of processes involved. An example was 
the use to constrain double beta decay rates [56|. The problems regarding approximations 
leading to effective local interaction forms of the optical potential are accentuated in all such 
analyses with the radial properties of the distorted waves from a nonlocal interaction being 
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quite different through the nuclear transition interaction region (interior and surface of the 
nucleus) from those of an "equivalent" local optical potential with which a good fit to elastic 
scattering data was found. 

While this study has centred upon the problems associated with an NA optical potential 
and its use in inelastic scattering analyses, many of the same issues attend analyses of 
nucleus-nucleus scattering data and conclusions about structure drawn from them. 
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